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We investigate the possibility that ultra-high energy cosmic rays (UHECRs) originate from the annihilation of 
relic superheavy (SH) dark matter in the Galactic halo. In order to fit the data on UHECRs, a cross section of 
(aAv) ~ 10-^^cm2(Mx/10^^GeV)^/^ is required if the SH dark matter follows a Navarro-Frenk- White (NFW) 
density profile. This would require extremely large-Z contributions to the annihilation cross section. An interesting 
finding of our calculation is that the annihilation in sub-galactic clumps of dark matter dominates over the 
annihilations in the smooth dark matter halo, thus implying much smaller values of the cross section needed to 
explain the observed fluxes of UHECRs. 



1. Introduction 

The detection of particles with energy in ex- 
cess of 10^° eV is one of the unsolved mysteries of 
modern astrophysics and might represent one of 
the many signals of the existence of physics well 
beyond the standard model of particle physics. 
In ||l|,^ it was proposed the appealing possibil- 
ity that this mystery might be related to another 
problem, the origin of dark matter. In fact the 
dark matter might be made of SH relics of the 
big bang with lifetimes possibly exceeding 

the age of the universe. The rare decays of these 
SH particles (with masses in excess of 10^^ GeV) 
naturally generates UHECRs from the top, that 
is as end-product of their decay. In order to satu- 
rate the dark matter content of the universe and 
at the same time explain the fluxes of UHECRs, 
the lifetime of the SH relics should be of order 
10^^ years, which is difficult to achieve for such 
massive particles unless some discrete conserved 



symmetry is introduced that prevents the natural 
decay of the relics on much shorter time scales. 

In this paper we investigate an alternative pos- 
sibility, consisting in the annihilation of SH relics 
rather than in their decay. The simple assump- 
tion that dark matter is a thermal relic limits the 
maximum mass of the dark matter particle to less 
than a few hundred TeV. However, it has been re- 
cently pointed out that dark particles might have 
never experienced local chemical equilibrium dur- 
ing the evolution of the universe, and that their 
mass may be in the range 10^^ to 10^^ GeV, much 
larger than the mass of thermal WIMPs 

In this paper we calculate the expected flux of 
UHECRs from the smooth dark matter content 
of our Galaxy, and the flux due to clumps of dark 
matter in the halo. 
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Figure 1. UHECR spectra from SH particle an- 
nihilation (upper panel) or decay (lower panel). 
For both figures the solid (dashed) lines are for 
Mx = 1012 GeV {Mx = lO^GeV). Cross sec- 
tions and decay rates are indicated in the plot. 



2. SH Dark matter annihilation and UHE- 
CRs 

In calculating the UHE cosmic ray flux from a 
smooth SH dark matter distribution in the galac- 
tic halo, we assume a halo density spherically 
symmetric about the galactic center, in the form 
of a NFW profile @ nxid) = ^^a+Xy' ' ""^^^ 
use ds = 3dQ = 24 kpc in our numerical esti- 
mates, where d© ~ 8 kpc is the distance of the 
solar system from the galactic center. The dimen- 
sionless parameter A^o may be found by requiring 
that the total mass of the Galaxy is 2 x lO^^M©. 
For simplicity, we will assume that each annihi- 
lation produces two quark initiated jets, each of 
energy Mx , while the decay of a SH particle pro- 
duces two jets, each of energy Mx/2. The en- 
ergy spectrum of observed UHE cosmic ray events 
from annihilation is determined by the energy 



distribution of the particles in each jet, which 
we take to equal the fragmentation function in 
the modified leading logarithmic approximation 
(MLLA) which was also employed in 

We also use the MLLA limiting spectrum in 
the results of Fig. |l|. Calculating the resulting 
UHE cosmic ray flux in the annihilation model, 
and comparing it to the similar calculation in the 
decay scenario, we obtain the results shown in 
Fig. H (the data points are from AG AS A). In or- 
der to provide enough events to explain the ob- 
served UHE cosmic rays, {aAv) has to be in the 
range 10"^^ cm^ to 10"^^ cm^, which is well in ex- 
cess of the unitarity bound to the Z-wave reaction 
cross section pO|jr^ , p^ . The unitarity bound es- 
sentially states that the annihilation cross section 
must be smaller than M^^. However, as empha- 
sized by Hui iQ , there are several ways to evade 
the bound. The annihilation cross section may 
be larger if there are fundamental length scales 
in the problem larger than M^^. 

A related issue is the typical energy of the an- 
nihilation products. In this paper we assume that 
annihilation produces two jets, each with energy 
approximately Mx- It is easy to imagine that 
with the finite-size effects discussed above, there 
is the possibility that annihilation will produce 
many soft particles, rather than essentially two 
particles each of energy AIx- An example that 
suits our needs is the annihilation of a monopole- 
antimonopole pair, as discussed at length in [Tl] ]. 
We regard the requisite size of the annihilation 
cross section to be the most unattractive feature 
of the annihilation scenario. 

This problem becomes less severe when the 
clumped component of the dark matter distribu- 
tion in the Galactic halo is taken into account. 
The spatial and mass distribution of the clumps 
is taken from numerical N-body simulations (see 
pT[ | and ||l^ ) . The presence of cuspy density pro- 
files in the clumps makes their contribution to the 
diffuse flux of UHECRs dominant over the NFW 
smooth component. The ratio of the clumped 
to smooth contributions for isothermal and NFW 
profiles for the clumps and for different values of 
Mx is plotted in Fig. 2. In the case of isother- 
mal clumps, where the density profile of a single 
clump scales like l/r^, the main contribution to 
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Figure 2. The ratio of events from the subclump 
component to the smooth component assuming 
either an isothermal or a NFW profile for the sub- 
clumps. 



the annihilation comes from the inner regions of 
the clumps, so that it is crucial to define a min- 
imum radius Rmim determined by the annihila- 
tion eSiciency. In |Tl| it is explained in detail the 
effect of two different definitions for Rmin- 

Another peculiar prediction of the annihila- 
tion of SH relics in clumps of dark matter in the 
halo consists of a peculiar pattern of small scale 
anisotropics. Several multiple events within de- 
gree scales would be observed if annihilation in 
the clumps generates the UHECRs. A detailed 
simulation was carried out in The model 

also predicts a large scale anisotropy in the direc- 
tion of the galactic center, which is an important 
test to be carried out by the future experiments, 
that will be able to see the direction of the galac- 
tic center. 

3. Conclusions 

Our calculations can be summarized in the fol- 
lowing points: 1) the annihilation scenario for 
the origin of UHECRs requires cross sections that 
may be larger that the unitarity bound. Although 
not a killing factor, for the reasons explained 
above, this is certainly an unappealing feature of 
the model; 2) the annihilations in clumps of galac- 



tic dark matter dominate over the contribution 
of a smooth NFW dark matter profile, requiring 
correspondingly smaller values for the annihila- 
tion cross section; 3) the pattern of anisotropics 
is such that both a large scale (dipole) anisotropy 
and a degree scale anisotropy is to be expected 
if UHECRs are generated in the annihilations of 
SH relics in galactic SH dark matter. 
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